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Abstract 



The transfer operator due to Bogomolny provides a convenient method for ob- 
taining a semiclassical approximation to the energy eigenvalues of a quantum 
system, no matter what the nature of the analogous classical system. In this 
paper, the method is applied to integrable systems which are rotationally in- 
variant, in two and three dimensions. In two dimensions, the transfer operator 
is expanded in a Fourier series in the angle variable, while in three dimensions 
it is expanded in spherical harmonics. In both cases, when the Fourier coef- 
ficients are evaluated using the stationary phase approximation, we arrive at 
the Einstein-Brillouin-Keller (EBK) quantization conditions. The associated 
Maslov indices are shown to agree with the results calculated by well-known 
simple rules. The theory is applied to several rotationally invariant systems, 
including the hydrogen atom and the isotropic harmonic oscillator in two and 
three dimensions, the circle billiard, a billiard inside a spherical cavity, and a 
harmonic potential with a singular magnetic flux line. 
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I. INTRODUCTION 



In furthering our understanding of the relationship between classical mechanics and 
quantum mechanics, semiclassical approximations play an important role. Periodic orbit 
theory, developed by Gutzwiller, Balian and Bloch, Berry and others employs the 

periodic orbits of the classical system to obtain a semiclassical approximation to the density 
of states or to individual energy eigenvalues of the analogous quantum system. Although 
formally elegant and satisfying, the theory is usually hard to apply in practice, because the 
periodic orbit sum is not absolutely convergent, and because it is difficult to find the periodic 
orbits in a systematic way. Alternative semiclassical approximations, which do not depend 
on knowing the periodic orbits, have been proposed by Bogomolny j|,[5| and by Doron 
and Smilansky |3[[7|. Exploiting the duality between the classical dynamics of a billiard 
inside a bounded region and the scattering of external particles by the system's boundary, 
Doron and Smilansky obtained semiclassical energy eigenvalues for the billiard system by 
constructing a semiclassical approximation to the scattering matrix for the exterior problem. 
In the theory proposed by Bogomolny, one chooses a Poincare surface of section (PSS) which 
is frequently crossed during the motion of the system, and one constructs a semiclassical 
transfer operator from the classical trajectories which take the system from one position on 
the PSS to another. For billiard systems, the approaches based on the scattering matrix 
and on the transfer operator can be shown to yield the same determinantal equation for the 
energy eigenvalues of the interior system ||||. 

The present paper had its origin in trying to derive the correct quantum energy eigen- 
values of the hydrogen atom (in three dimensions) by means of Bogomolny's semiclassical 
transfer operator. By adding a small 1/r 2 term to the Coulomb potential — a device which 
prevents the transfer operator from being singular — we achieved this goal. However, an 
essential part of the derivation was the use of the stationary phase approximation. Sub- 
sequently, we generalized our approach and showed that it led to the well-known Einstein- 
Brillouin-Keller (EBK) quantization rules PHT?1- When viewed in this light, there is nothing 
special about our solution for the hydrogen atom. In fact, the main result of our paper can 
be summed up concisely: for rotationally invariant integrable systems, Bogomolny's transfer 
operator, plus the stationary phase approximation, yields EBK quantization. In a separate 
publication JTB] we plan to show how the EBK quantization rules can also be obtained from 
Bogomolny's transfer operator formulated in terms of the angle-action variables. 

The plan of the paper is as follows. After a brief description of the transfer operator 
in the next section, in Sec. f ll| we construct the transfer operator for a two-dimensional 
system having circular symmetry. By making a suitable Fourier expansion and evaluating 
the Fourier coefficients by means of the stationary phase approximation, we derive the EBK 
quantization conditions. This general formulation is applied in Sec. [IV] to the hydrogen atom 
(plus 1/r 2 potential), the circular harmonic oscillator (plus 1/r 2 potential), a harmonic 
potential plus a singular magnetic flux line, the circle billiard, and the annulus billiard. 
In Sec. [V] a similar approach is described for three-dimensional systems having spherical 
symmetry. The resulting EBK quantization conditions are applied in Sec. [VT| to the hydrogen 
atom (plus 1/r 2 potential), the isotropic harmonic oscillator (plus 1/r 2 potential), and a 
billiard inside a spherical cavity. The paper concludes with a discussion of our results for 
the EBK energy eigenvalues in comparison with the exact quantum energies. 
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II. THE TRANSFER OPERATOR AND THE DETERMINANTAL EQUATION 



We begin with a brief description of Bogomolny's semiclassical transfer operator. For 
a system with / freedoms, the PSS in configuration space is a surface or hypersurface of 
dimension / — 1, and the transfer operator in the coordinate representation is H 



T(q",q';E) 



y I 

k. (2*1/0 C/-U/a 



det 



d 2 S(q", q'; E) 



dq"dq' 



exp[iS{q",q';E)/h-ifnr/2}, (2.1) 



where q' and q" denote / — 1 generalized coordinates for two points located on the PSS. The 
summation is over all classical trajectories which go from q' to q", crossing the PSS at these 
points in the same sense and at no other points (in the same sense) in going from q' to q" . 
For each such trajectory one needs the action at energy E, denoted by S(q", q'\ E), and the 
phase index /x, which is related to the occurrence of caustics — points on the trajectory at 
which the semiclassical approximation is not valid. The matrix of second derivatives of the 
action has dimension / — 1. 

If T(E) is the transfer operator and I is the unit operator, the corresponding semiclassical 
energy eigenvalues of the quantum system are determined from the condition M 



det[I-T(E)) = 0. 



(2.2) 



When properly formulated, the T-operator is unitary ||. In the past few years there have 
been several applications of the transfer operator based on constructing an approximation 
to T(E) in coordinate space H-El. However, when the system being treated has rotational 
symmetry, it is better to treat the transfer operator in the angular momentum representation, 



since it is then diagonal |]I9 |. In the next section we show how this may be carried out for 
two-dimensional systems with circular symmetry. 



III. TWO-DIMENSIONAL SYSTEMS WITH CIRCULAR SYMMETRY 

Let us consider a particle of unit mass moving in two dimensions in a potential V(r). 
The Hamiltonian is 

H = P j + §- 2 +V(r) (3.1) 

where p r = r and p^ = r 2 <fi are the momenta conjugate to the polar coordinates r and 
4> describing the particle's position. Since the angular momentum p^ is a constant of the 
motion, let us denote it as L. For given E and L, the turning points of the classical motion 
along the radial direction are determined by 

^T 2 +V{r) = E. (3.2) 

Clearly, the turning point radii, r_ and r+, depend on L as well as E. 

In setting up the transfer operator, we choose the Poincare surface of section (PSS) to be 
a circle of radius R. While in principle any radius between r_ and r + could be used, there 
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is a natural choice for given energy E. At that energy there is a trajectory with maximal 
angular momentum L max (E) which is a circle, corresponding to the radial kinetic energy 
being zero. We define R to be the radius of this circle. It can also be thought of as the 
circle for which the radial turning points coincide. Then, trajectories at energy E having 
\L\ < \L max (E)\ have nonzero radial kinetic energy, and, therefore, must repeatedly cross 
this circle, making it a suitable choice for the PSS. (It is easy to show from equation ( |3.2| ) 
and the derivative of this equation that R is the solution of 2V(r) + rV'(r) = 2E.) 

The coordinate q on the PSS will be taken to be the polar coordinate 0. From equation 
( ^.lp the transfer operator from to 0' (where < < 2n and < 0' < 2n) on the PSS is 

1/2 



(2vr^) 1 /2 



«9 2 ^(0',0;£) 



d(f)'d(J) 



exp[iSj (<f)', 4>;E)/h — ifijir/2], (3.3) 



where j labels different possible classical trajectories at energy E which go from to 0' 
without crossing the PSS (in the same sense) at any other point. Let us define the angle 7 
to be 7 = 0' — (modulo 2tt) in order that < 7 < 2tt. The trajectories can be labelled in 
such a way that the angle traversed by the particle in going from to 0' is = 7 4- 2irj, 
where < 7 < 2tt and j is an integer (positive, negative or zero). The action along the jth 
trajectory can be written as S(£W; E). Also, we can write 

d 2 S(^;E) _ d 2 S(^;E) 

9090' <97 2 1 ' 1 

Thus, in accord with the invariance of the Hamiltonian under rotations about the origin, the 
transfer operator can be expressed in terms of the relative angle 7 only. We can, therefore, 
expand (|3.3|) in a Fourier series: 



T(0',0;£) = T( 7 ;£)= ]T C m (E) exp(*m 7 ), (3.5) 

m=— 00 

where the expansion coefficients are 

1 r 27T 

C m {E) = — T( T ,E)exp(-im 1 )d 1 . (3.6) 

Z7T JO 

We now construct a matrix representation of the transfer operator using the basis 
{(27r)~ 1//2 exp(im0)}. Since 7 = 0' — (modulo 2tt), a typical matrix element is 

T mim2 (E) = — / d(f) d0'exp(-zmi0')T(7;£;)exp(im20) = 2n C mi (E)5 mim2 . (3.7) 
2w Jo Jo 

Thus, the T-matrix is diagonal in this representation, and its eigenvalues (as a function of 
E) are just the diagonal elements. (Note that these T-matrix eigenvalues should not be 
confused with the semiclassical energy eigenvalues.) Denoting the mth eigenvalue curve as 
X m (E), we obtain from equations ( |3.3|) to 



r2ir 

X m (E) = / T(7; E) exp(—im'f)d'y 







(2mh) 



1 f 2lT 
j— ^exp(-z/i i 7r/2) ^ d 7 



d 2 S(^;E)\ 1/2 



<9 7 2 



exp[iS(£ ij) ;E)/h-im-f}. (3.8) 
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We remind the reader that ft) depends on 7 through the definition ft) = 7 4- 2nj. 

Up to this point we have made no approximations other than the approximation involved 
in deriving Bogomolny's semiclassical transfer operator. We now evaluate the integrals in 
Q3.8| ) using the stationary phase approximation. For the jth integral the point at which the 
phase is stationary is determined by the equation 

8g (f^' =rofi . (3.9) 

The left-hand side of this relation is the classical angular momentum L for the jth trajectory. 
Thus, the stationary phase condition effectively quantizes the angular momentum of the 
particle. Because the possible trajectories from to 0' at energy E are uniquely specified by 
the angular momentum, it is clear that for each value of m (positive, negative or zero) there 
is at most one trajectory satisfying ( [3.91) . (A solution exists if and only if \L\ < \L max (E)\.) 
The value of j for this trajectory will be denoted j m . Thus, for given m, we denote the 
solution of equation (|3.9| ) (when it exists) as 7 m , and the corresponding angle traversed by 
the particle in going from to 0' as £ m = 7 m + 2irj m . 

When a solution 7 m of equation Q3.9|) exists for a particular value of j m , we can evaluate 
the integral in the usual way, assuming that d 2 S/d r y 2 is a relatively slowly varying function 
of 7. Introducing the symbol v m through the definition 

v m = if 




u m = l if l w ) <0 (3.10) 

\ I / 7=7 m 

and henceforth denoting //,■ for the trajectory j m as /x m , we obtain 

\ m {E) w exp[iS(£ m ; E)/% - im^ m - i([i m + v m )ir/2]. (3.11) 

The fact that these approximate eigenvalues of the T-matrix have unit modulus is consistent 
with the T-matrix being unitary [|J. Note that the symbols \i m and v m have the same 
meaning as in the paper by Creagh, Robbins and Littlejohn |24| . 



It is useful to further simplify this expression by splitting the action for the trajectory 
j m into radial and angular parts. The angular part, evaluated at j m , may be denoted by 

SangiZm', E) = J p^dcj) = mU{^ m + 27TJ m ), (3.12) 

and the radial part, evaluated at the angular momentum L = mh determined by the sta- 
tionary phase condition, is 

/f r + 
\p r \ \dr\ = 2 / \pr\dr. (3.13) 

Using the fact that exp(im2irj m ) = 1, we arrive at the expression 

X m (E) ps exp[iS rad (L = mh; E)/% - z(/i m + v m )-n/2}. (3.14) 
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The semiclassical energy eigenvalues of the quantum system are found from the deter- 
minantal equation ( |2.2|) , which is satisfied whenever an eigenvalue of the T-matrix is equal 
to unity. Thus, the condition for an energy eigenvalue is that \ m (E) = exp(i2im r ) . From 



equation (|3.14j ) this yields 



S T ad(L = mh; E) = 2iih(n r + cr m /4) n r = 0, 1, 2, . . . (3.15) 

where we define the Maslov index a m = fi m + v m . This is associated with a complete cycle 
of the radial motion for a trajectory at energy E and angular momentum L = mh. The 
allowed values of n r in ( |3.15| ) are determined by the assumption that SVad > 0. 



In Appendix |A] it is shown that, for smooth potentials, the combination o m = fi m + v m 
is always equal to 2. This is the result one would obtain for the Maslov index in EBK 
quantization |TT| , |T2"11 using the simple rule of counting 1 for each of the soft turnarounds 
during a complete cycle of the radial motion. It is also shown in Appendix [A] that if the 
particle is confined inside a circular disk with a hard wall (Dirichlet boundary condition on 
the wave function), the result for a m is 3. This agrees with the simple rule for computing 
the EBK Maslov index by counting 1 for the soft turnaround at the inner radial turning 
point and 2 for the collision with the disk boundary. Thus, for systems having circular 
symmetry in two dimensions, Bogomolny's transfer operator (modified using the stationary 
phase approximation) leads to EBK quantization, with the Maslov index for the radial 
motion computed by the well-known simple rules. Note that canonical invariant, 

even though the transfer operator is not canonically invariant (since it depends on the choice 



of the PSS). In addition to ( |3.15| ), the other EBK quantization condition, L = mh, was, 
of course, obtained from the stationary phase condition ( |3.9| ). It is worth noting that for 
systems that are invariant under time reversal, the energy eigenvalues with m^O are doubly 
degenerate. 



IV. APPLICATION TO SYSTEMS IN TWO DIMENSIONS 
A. The Coulomb plus 1/r 2 potential 

Let us now apply this general formulation to the hydrogen atom in two dimensions. In 
fact we shall treat a slightly more complicated potential, namely the Coulomb potential 
plus a term proportional to 1/r 2 . For the pure Coulomb potential, there is a one-parameter 
family of ellipses which start out from a given point on the PSS and return to the same 
point. Thus, this point is a focal point. By adding the 1/r 2 term, we ensure that the 
trajectories are not ellipses, and thereby avoid difficulties associated with the initial point 
being a conjugate point. 

Assuming that the nucleus is stationary at the origin, we take the potential to be 

V ^ = ~\ ± ^- (4-1) 

For convenience we have taken the electronic charge e to be unity, and we have written the 
strength of the 1/r 2 potential as a 2 /2, where a has the dimensions of angular momentum. 
The 1/r 2 potential may be repulsive or attractive, and provided a 2 is not too large in the 
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repulsive case, the electron will always be bound to the nucleus, implying that the energy 
E is negative. From equation ( |3.2|) , the classical turning points of the radial motion occur 
at radii r_ and r + given by 



(3=[l-2\E\(L 2 ±a 2 )] 1 / 2 . (4.3) 

The radius of the Poincare circle defined in the previous section is R = 1/(2\E\) (although 
we shall not make explicit use of this in what follows). For the repulsive 1/r 2 potential, each 
trajectory traverses an angle less than 2ir before returning to the PSS, while in the attractive 
case, the trajectories go through angles greater than 2tt. This means that the possible 
trajectories are qualitatively different in the two cases, there being only two possibilities in 
the former case and many in the latter case. Note from (|4.3| ) that, for given E, the maximum 
possible value of the classical angular momentum is given by 

[W£)] 2 = ^T« 2 . (4.4) 

To obtain the semiclassical energy eigenvalues from equation (|3.15|) , we must calculate 
the radial action integral with L = mh, as in (|3.13|) . Using equation (|3.1|) with L = mh to 
solve for \p r \ as a function of r, we obtain 

7 / 2 \ 1/2 

S iad (L = mh; E) = 2 J \ Pr \dr = vr ^— j - 2n(m 2 h 2 ± a 2 ) 1 ' 2 (4.5) 

Here the first term on the right-hand side is the action of each member of the family of 
elliptical orbits at energy E of the pure Coulomb potential. To write down the EBK quanti- 
zation condition from equation fl3.15| ), we set a m = 2, corresponding to two soft turnarounds 
at the radial turning points (see Appendix |A|). Hence, 



/ 2 \ 1 

-\—\ - 2n(m 2 ± a 2 /h 2 ) 1/2 = {2n T + 1)tt, n P = 0,1,2,... (4.6) 



Using the fact that E is negative for the bound state solutions we are considering, we obtain 

1 



E, 



""" ' 2h 2 [n r + \ + (m 2 ± a 2 /h 2 y/ 2 } 2 



2 

m = 0,±l,±2, ... n r = 0,1,2,... (4.7) 

This expression gives the approximate semiclassical energy eigenvalues for the Coulomb 
plus 1/r 2 potential. Note that the energy is the same for positive and negative values of m, 
implying that the energy eigenvalues are doubly degenerate for m ^ and nondegenerate 
for m = 0. The allowed values of m are constrained by the condition \m\h < \L max (E)\, 
with \L max (E)\ given by equation 
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The pure Coulomb potential is obtained by letting a — > in equation ( |4.7|) . Until now 
we have assumed that a is nonzero and sufficiently large (presumably, a ^> Ti) that the 
trajectories are not close to being ellipses, thereby avoiding the point 0' on the PSS being 
a focal point. At this stage, however, it is permissible to relax this requirement and let a 
become zero. Putting n = \m\ + n r , we obtain 



1 



E n = - n*2t„ , ly, . n= 0,1,2,... (4.8) 



(It is satisfying that the limit does not depend on the sign of the a 2 /(2r 2 ) term in the 
potential, despite the fact that the two problems are quite different, as mentioned earlier.) 
This expression for the energy eigenvalues is exactly the same as the result found by solving 
the two-dimensional Schrodinger equation for the Coulomb potential. The eigenvalues do not 
depend explicitly on m, but it is clear from the definition of n that \m\ < n. This condition, 
which also arises, for example, in solving the radial Schrodinger equation by the method 
of series expansion, or by using group-theoretical considerations, correctly determines the 
degeneracies of the energy levels given by ( |4.8| ) . 

In an earlier study of the hydrogen atom in two dimensions using Bogomolny's transfer 
operator [J22|J, the PSS was chosen to be a radial line. The "half-mapping" transfer operators 



introduced by Haggerty [[18] were used to avoid the problem associated with the family of 
trajectories (ellipses) at energy E starting from a point on the PSS and returning to the 
same point on the PSS. The outcome of this work was similar to equation (|4.8|) but with 
(n+ \) 2 replaced by (n+ |) 2 . We now realize that this peculiar result was due to an incorrect 
assignment of the phase indices associated with the elliptical trajectories. In Fig. 2 of Ref. 



22] , there is a caustic associated with the longer solid trajectory, but there is no caustic 



associated with the shorter solid trajectory. When this fact is properly taken into account, 



the energy levels turn out to be the same as equation (|4.8| ). Thus, the energy levels of this 



system do not depend on the choice of the PSS, at least for the two choices considered. 



B. The circular harmonic oscillator plus 1/r 2 potential 
A particle moving in two dimensions in a circular harmonic oscillator plus 1/r 2 potential 



has been treated in an earlier paper [23| using a slightly different method based on Bo- 



gomolny's transfer operator. Here we show that our general formulation of Sec. |TJ leads 
quickly to the same results for the energy eigenvalues of this system. 
We take the potential to be 

1 a 2 

V(r) = -u J 2 r 2 ±—, (4.9) 

where, as in the last subsection, ct 2 /2 is the strength of the 1/r 2 potential, which may be 
attractive or repulsive. When the particle has energy E, equation (|3.2|) leads to the following 
expression for the classical turning-point radii: 

rl = £± ' E2 '^f ±a2 >]' /2 . (4.10) 



s 



Using this, and solving equation to find \p r \ as a function of r, one finds, 

S rad (L; E) = 2 r \ Pr \dr = — - tt(L 2 ± a 2 ) 1 ' 2 . (4.11) 



When this is substituted in equation (|3.15|) , with L set equal to mh and a m set equal to 
2 (corresponding to soft turnarounds at r_ and r+), we obtain for the energy eigenvalues 
belonging to a given value of m 

E mnr = hu[2n r + (m 2 ± a 2 /h 2 f/ 2 + 1], 
m = 0, ±1,±2, ... n r = 0,1,2... (4.12) 

This expression agrees with the result obtained from an exact solution of the Schrodinger 
equation (see, for example, Flugge IffiD- A plausible explanation of this agreement, despite 
the approximations inherent in our semiclassical approach, has been put forward earlier 
2~3fl and will be briefly discussed in the final section of the paper. Note that the energy 



eigenvalues are doubly degenerate when m^fl and nondegenerate when m = 0. 



C. Circular harmonic oscillator plus singular magnetic flux line 

The development in the last subsection can be extended to include a singular magnetic 
flux line passing through the origin. This means that, in addition to the circular harmonic 
oscillator potential and the 1/r 2 potential, the particle motion occurs in the presence of a 
magnetic field (perpendicular to the plane of the motion) having the form of a (^-function 
singularity at the origin. The magnetic field breaks the time-reversal symmetry and removes 
the degeneracy of the energy eigenvalues when m^fl. 

It is convenient to paramaterize the strength of the flux line by the positive quantity 
5 = e$/(hc), where e is the magnitude of the charge on the particle and $ is the total 
magnetic flux through the singular point. In the presence of the flux line the Hamiltonian 
is (see Brack et al. p6| ) 

H = i + <n^? +v(r) (413) 

where p^, the momentum canonical to the coordinate 0, is a constant of the motion, and 
V(r) is given by equation ( |4.9| ). Then, for the trajectory labelled by j in equation (|3.8|), the 



angular part of the action is 

S ang (£ (j) ; E) = P< £® = p^ + 2nj). (4.14) 



Thus, the stationary phase condition ( |3.9| ) becomes p^ = mh. However, by comparing ( 4.13|) 
with equation (3.1) we see that (p^ — 5) 2 now plays the role that L 2 played in Sec. |TTT|. For 



given m the quantization condition ( |3.15| ) becomes 

S Iad (L= \mh-5\;E) = 2irh(n r + cr m /4), n r = 0, 1, 2, . . . (4.15) 



Making this replacement in ( |4.11| ) and setting a m = 2 (corresponding to soft turnarounds 



at r__ and r + ), we obtain the following expression for the energy eigenvalues: 
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E mnr = hu{2n T + [(m - 6/h) 2 ± a 2 /^] 1 / 2 + 1}, 

m = 0,±1,±2,... n r = 0,1,2,... (4.16) 

For the case of a singular flux line at the centre of a circular harmonic oscillator potential, 
we can let a — > in ( f4.16| ). The resulting energy eigenvalues are 

Emn r = huj(2n r + \m — 5/h\ + 1), 
m = 0, ±1,±2, ... n r = 0,1,2,... (4.17) 

which clearly are different for positive and negative values of m. This result is the same 
as the exact analytical solution of the Schrodinger equation for this problem, which can be 
obtained from the exact harmonic oscillator solution by replacing \m\H by \m% — 8\ wherever 
it occurs p6 . 



D. The circle billiard 

In this section we apply the general formulation based on the transfer operator to a 
particle moving in a constant potential (which we take to be zero) inside a circle of radius 
R. We show that this leads to well-known results for the EBK energy eigenvalues. 

At energy E and angular momentum L, the radial motion has an inner turning point at 
the radius r_ given by \L\ = r_(2S) 1//2 . Choosing the radius of the Poincare circle to be R 
(or just slightly less than R so that the PSS is crossed just after the particle has made a 
collision with the boundary) ensures that all trajectories cross the PSS. 

At fixed values of E and L, the radial part of the action integral is 

SUL; E) = 2 J r R \p r \dr = 2 J r R (2E - L 2 /r 2 f' 2 dr 
= 2{2ER 2 - L 2 ) 1 ' 2 - 2\L\ cos^ 1 [|L|/(2 J E) 1 / 2 J R]. (4.18) 

Writing E = h 2 k 2 /2 and setting L = mK, with the values of m restricted by the condition 
\iTi\h < \L max (E)\ = (2Ey/ 2 R, we obtain from equation (|3.15 ) the following condition for 
an approximate energy eigenvalue of the quantum system: 



{k 2 R 2 - m 2 fl 2 - \m\ cos- 1 [|m|/(fci?)] = 7i{n r + 



-) 



m = 0,±l,±2,... n r = 0,1,2,... (4.19) 

Here we have put a m = 3, corresponding to a soft turnaround at the inner turning point 
and a hard- wall collision at the circle boundary. (See the discussion in Appendix [A].) Equa- 
tion ( |4.19| ) can be solved numerically to determine the EBK energy eigenvalues. (This is 
equivalent to finding the zeros of the Bessel function J m (kR) when it is approximated by 
the leading term of the Debye asymptotic expansion. See Ref. 0, p 336.) Results for the 
lowest energy eigenvalues have been tabulated by Keller and Rubinow JEU and by Brack 
and Bhaduri (see Ref. ||, p 88). The fractional difference between the EBK eigenvalues and 
the exact energy eigenvalues (determined from the zeros of the Bessel function J m (kR)) was 
found to decrease fairly rapidly with increasing energy. 
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E. The annulus billiard 



The annulus billiard consists of a particle moving in a constant potential (which we take 
to be zero) in the region between two concentric circles of radii R and a. (We assume 
R > a.) In this section we show how equation ( [3.15| ) may be used to obtain the EBK energy 
eigenvalues for this system. 

As in the case of the circle billiard, we choose the Poincare circle to have radius R. At 
energy E and angular momentum L the radial motion may have an inner turning point at 
the radius r_ given by \L\ = r_(2E) 1 / 2 . Provided a < r_, the minimum value of r during 
the radial motion will be r_. However, if a > r_, the radial motion is reversed by a hard- wall 
collision at r = a. In the former case with a soft turnaround at r_, the Maslov index is 
a m = 3, as for the circle billiard. In the latter case, there are two hard-wall collisions, and 
the Maslov index is a m = 4. 

The radial part of the action involves the same integral as in equation Qllfc ), but 



is now evaluated at the limits r min and R, where r min is r_ or a, whichever is larger. 
Putting E = h 2 k 2 /2 and setting L = mh, with the values of m restricted by the condition 
\m\% < \L max (E)\ = (2E) l l 2 R, we obtain the following condition for an approximate energy 
eigenvalue of the quantum system, similar to equation (|4.19|) : 



(k 2 r 2 - m 2 ) 1/2 !L„ - H cos-^lmlAHllLn = ^ + <W 4 )> 

m = 0, ±1,±2, ... n r = 0,1,2,... (4.20) 

Here, a m = 3 if r min = r_, and a m = 4 if r min = a. The solution of this equation 
gives the EBK energy eigenvalue E mrLr . Numerical values for the lowest 30 distinct EBK 



energy eigenvalues are tabulated in Ref. [28] for three cases: a = 0.1R, a = 0.3R, and 



a = 0.5R. Also tabulated in Ref. |28[] are the energy eigenvalues computed from Bogomolny's 
transfer operator without making use of the stationary phase approximation or any other 
approximation. The results differ appreciably from the EBK energy eigenvalues. This 
draws attention to the fact that the stationary phase approximation used to derive equation 
( P-15| ) constitutes an approximation in addition to the main semiclassical approximation 
contained in Bogomolny's transfer operator. It is clear that Bogomolny's approach and 
EBK quantization are generally equivalent only to leading order in %. 



V. THREE-DIMENSIONAL SYSTEMS WITH SPHERICAL SYMMETRY 

A theory of the transfer operator for three-dimensional systems with spherical symmetry 
is only slightly more complicated than in two dimensions. For a particle of unit mass moving 
in a potential V(r), the Hamiltonian is 

where p r = r, pg = r 2 9 and = r 2 sin 2 9(j) are the momenta conjugate to the spherical 
polar coordinates (r, 6, <fi) describing the particle's position. Because the component of the 
angular momentum about the polar axis, p^, is a constant of the motion, let us denote it 
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as L z . In addition, the square of the total angular momentum, L 2 = 
constant of the motion. Therefore, the Hamiltonian may be written as 

H = — H + V(r) 

2 2r 2 v ; 

For given and L, the turning points of the radial motion are determined by 



+ L\j sin 2 9, is 



2r< 



+ V(r) = E. 



(5.2) 



(5.3) 



Thus, as in two dimensions, the turning point radii, r_ and r + , depend on L 2 as well as E. 

In constructing the transfer operator, we choose the PSS to be a sphere of radius R. The 
argument following equation ( |3.2|) applies equally well to the three-dimensional situation. 
Therefore, we choose R to be equal to the radius of the circular trajectory at energy E 
for which the angular momentum has its maximum possible value. Then, for any angular 
momentum \L\ < L max (E), the trajectories of the system will repeatedly cross this surface. 

The generalized coordinates q for the transfer operator are the polar angles (8, <p) on the 
Poincare sphere, which we denote collectively as Q. Then, from equation ( |2.1| ) with / = 3, 
the transfer operator from Q to Q' is 



T(n',n;E) = J2 



cl.tr. 



2irih 



det 



d 2 S 



dttdtt* 



exp[iS(fl', Q; E)/h — ifnr/2], 



(5.4) 



Here the sum is over all possible classical trajectories at energy E which go from Q to 
Q' without crossing the PSS (in the same sense) at any other point. Let 7 be the angle 
subtended at the origin by these points, with < 7 < tt. All possible trajectories from Q to 
Q' lie in the plane defined by these points and the origin, and each such trajectory can be 
uniquely identified by the total angle traversed as the particle moves along the trajectory. 
The trajectories fall into two classes, which we treat separately. In the first class, the angles 
traversed are = 2irj + 7, j — 0, 1,2,.. ., while in the second class, the angles traversed 
are £_ = 2-kj — 7, j = 1, 2, . . .. The sum in ( [5.4|) is over the trajectories of both classes. 

The determinant of the second derivatives of S(£± , E) is evaluated in Appendix ||, where 
it is found that 



det 



d 2 S 



dVtdVt' 



sin 7 



dS 



d 2 S 



<9 7 2 



(5.5) 



Substituting in equation (|5.4|) , we obtain 



cl.tr. 



2-KlTl (SUI7) 1 / 2 



dS 



97 



d 2 S 



<9 7 2 



exp[zS(d J) ;£)A-^V2] 



OX 



(5.6) 



The spherical symmetry of the system implies that T(Q', Q; E) depends only on the angle 
7 subtended at the origin by the points Q and Q' on the Poincare sphere. Therefore, we can 
expand the transfer operator in a Legendre series: 



T(Q', Q- E) = T( 7 ; E) = Y, Q(E) P z (cos T ), 



(5.7) 
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where the expansion coefficients are 



cm 



2/ + i r 



rn 

/ TOy; E) Pi(cosy) sin 7 dy. 
Jo 



(5.8) 



The spherical harmonic addition theorem 



47T 



(5.9) 



m=—l 



allows us to write 



An 



T(tt, O; E) = T( 7 ; E) = J2 Ci(E)—j E ^(O^tf 



(5.10) 



where Yj m (fl) is a spherical harmonic function. We can now calculate matrix elements of 
the T-operator in the angular momentum representation. From the orthonormality of the 
spherical harmonics, a general matrix element is found to be 



T hmi>l2m2 {E) = fdnfdMY£ mi {n)T(M,n;E) Y hm2 (n') 

An 



2/i + 1 



Ci 1 (E)8i 1 i 2 8. 



m\m,2 • 



(5.11) 



Thus, the T-matrix is diagonal in this representation, and its eigenvalues (as a function of 
E) are just the diagonal elements. Denoting the /m-eigenvalue curve, which is (21 + l)-fold 
degenerate, as Xi m (E), we obtain from equations (p.6|), ( |5.8[) and ( |5.11| ), 



^lm(E) 



ih 



cl.tr. 



E L d ^ 



dS 



dy 



d 2 S 



dy 2 



exp[iS{^£; E)/% - P i (cos7)(sin 7 ) 1 / 2 . (5.12) 



.(J) 



1/2 



So far we have made no approximations beyond those used to derive Bogomolny's semi- 
classical transfer operator. For a given system, one can calculate the action S(^;E) for 
each possible trajectory, as well as the derivatives with respect to 7, and hence evaluate the 
integrals. However, to proceed in a manner analogous to our treatment of two-dimensional 
systems, we make use of the following asymptotic expansion for P;(cos7), valid for large 
values of I (see Ref. p7j , equations 8.10.7 and 6.1.46): 

. 1/2 



Pi (cos 7) 



1 



1 



(l + l) 1 ' 2 l27rsin7 



{exp[i(l + |)7 - i7r/4] + exp[-i(Z + \)y + in/A}}, 

(5.13) 



where terms of higher order in l/l have been neglected. The same approximation has been 
used in the problem of scattering of waves by a sphere to show how the limit of geometrical 
optics can be obtained from physical optics PD|. After inserting ( |5.13| ) in equation ( p.!2| ), 
we can evaluate the resulting integrals using the stationary phase approximation. 
The two terms in ( 5.13Q yield the stationary phase conditions 
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E) 



± {i + \)h = o, 



(5.14) 



where the ± between the terms refers to the two terms in the asymptotic expansion fl5.13| ). 
At this point we find it convenient to treat separately the trajectories that traverse angles 
= 2rtj + 7, j — 0, 1, 2, . . . and those that traverse angles £_ = 2-n-j — 7, j — 1, 2, . . ., 
where in both cases < 7 < n. We shall show that both classes of trajectory lead to the 
same stationary phase condition, but that for a given value of I, the corresponding trajectory 
is either in one class or the other. 

Let us suppose for the moment that the trajectory under consideration belongs to the first 



class. Because the action increases with the traversed angle £+\ we have dS/ 'd^=dS/ 'd£ 



> 



0. Then equation Q5.14 ) can be satisfied by taking the second term in ( |5.13| ). Since dS/dj 
is the magnitude of the total angular momentum, the stationary phase condition (|5.14[) 
effectively quantizes the total angular momentum: 



z = 0,1,2, 



(5.15) 



Furthermore, specifying the total angular momentum (with L < L max (E); see below equa- 
tion (|5.3|)) completely determines the trajectory at energy E in the plane containing Q, Q' 



and the origin. Thus, for a given value of I in ( 5.15 ), there is at most one trajectory, labelled 
by ji > 0, which contributes to the sum in equation (|5.12 ). For this trajectory (which, 
by assumption, belongs to the first class), we denote the angle traversed by the particle as 
6 — li + where (in the range between and 7r) is the stationary point determined by 



equation ( 5.14 ). We also denote the corresponding action by Si and the phase index by 
When slowly varying quantities are taken outside the integral (to be evaluated at 7 = 7^), 



the T-matrix eigenvalue curve ( |5.12j ) becomes 



OS 
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7=7; 



d 2 S 



<9 7 2 



7=7; 



x / dj exp[iSi(r, E) /h - i(l + ^7 - wr/2 + i7r/4] 



(5.16) 



The integral is now evaluated in the usual way. Introducing the phase index v\ as in 
equation ( |3.1(J| ), we obtain 



X lm (E) w exp[iSifa E)/h - i(l + |)7 - z(/i, + i/ { )tt/2] 



(5.17) 



As in the two-dimensional case, these approximate T-matrix eigenvalue curves have unit 
modulus, consistent with the T-matrix being unitary. 

We now split the action E) into radial and angular parts. From the stationary phase 
condition, equation ( [5. 14]) , the magnitude of the angular momentum along the trajectory is 
(I + h)h. Thus, since the particle traverses the angle ^1 — ^1 + 2irji when it moves along this 
trajectory, we obtain for the angular part of the action 



S^{i l] E)/h = {l + \){ ll + 27i Jl ). 



(5.18) 



The radial part of the action, evaluated for L 2 = (I + h) 2 h 2 , is denoted as 
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S rAd {L 2 = {l + \fh 2 ] E)= <f \p r \\dr\ = 2 £ + \p r \dr. (5.19) 

Hence, equation ( |5.17|) becomes 

X lm (E) « exp[^ rad (L 2 = (l + \fh 2 - E)/h - + 14 - 2^/2]. (5.20) 

Next we suppose that the trajectory corresponding to a given value of I belongs to the 
second class, in which case the angle traversed by the particle has the form = 2irj — 7, 
j = 1,2,... Here too the action must increase with the traversed angle £_ , implying that 
dS/d'y = —dS/d£_ < 0. Then the stationary phase condition, equation ( 5.14 ), can be 



satisfied by taking the first term in the asymptotic expansion ( [5.131 ). One obtains in this 
case the same condition, equation ( |5.15| ), for the magnitude of the angular momentum. 
However, the T-matrix eigenvalue curve corresponding to equation ( |5.17| ) becomes in this 
case 

X lm (E) « exp[iS,(&; E)/n + i(l + |)7i - i{m + ui + 1)tt/2], (5.21) 
where £j = 2ixji — 7; and ji > 1. The angular part of the action in this case is 

E )/ h = (* + §) (27TJ, - T{ ), (5.22) 

which yields the T-matrix eigenvalue curve 

Xim(E) « exp[^ rad (L 2 = (Z + ±) 2 fr 2 ; - + z/ z - 2j ; + 1)tt/2]. (5.23) 

The semiclassical energy eigenvalues of the quantum system are found from the deter- 
minantal equation fl2.2|) , which is satisfied whenever an eigenvalue of the T-matrix is equal 
to unity. Thus, the condition for an energy eigenvalue is that Xi m (E) = exp(z27rn r ). ^From 
equations (|5.20|) and ( [5.23 ) this yields the following condition for an energy eigenvalue: 



S rad (L 2 = il + \) 2 Ti 2 -E) = 27rh(n r + <r l /4), n r = 0,1,2... (5.24) 

where the Maslov index ai associated with a complete cycle of the radial motion at energy 
E and L 2 = {I + l/2) 2 h 2 is defined to be 

ai = fii + vi- 2ji for trajectories in class & = 2-rrji + j h j t = 0, 1, . . . 

&i — l^i + V\ — 2ji + 1 for trajectories in class & = 2nji — ji, ji = 1,2,... (5.25) 

As in the two-dimensional case, the values of n r in (|5.24j) are determined by the assumption 

that ^rad > 0. 

It is shown in Appendix [A| that this definition always leads to the result ai = 2 for 
soft potentials, and to o\ = 3 for a particle inside a spherical cavity, making a hard-wall 
collision with the boundary in each cycle of the radial motion. These results are consistent 
with computing o\ by the simple EBK rules of counting 1 for each soft turnaround of the 
radial motion and 2 for each hard-wall collision. Thus, equations ( 5.24|) and ( |5.25| ) give 



the EBK quantization condition for the radial part of the action in three dimensions. The 
other EBK quantization conditions are L 2 = (I + \) 2 Ti and L z = mh with —l<m<l. 
The first of these was obtained earlier from the stationary phase condition. The second 
is a natural interpretation of the quantum number m introduced through the expansion 
of the T-operator in spherical harmonics. Because of the spherical symmetry, the energy 
eigenvalues corresponding to a particular value of I are (21 + l)-fold degenerate. 
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VI. APPLICATION TO SYSTEMS WITH SPHERICAL SYMMETRY 



A. The Coulomb plus 1/r 2 potential in three dimensions 

As in the two-dimensional case, the pure Coulomb potential has the special property 
that all classical trajectories (ellipses) starting from a given point on the PSS will return to 
the same point, which is, therefore, a focal point. To avoid this singular behaviour we shall 
again add a small 1/r 2 term to the potential, which may be attractive or repulsive. Thus, 
we take the Hamiltonian corresponding to equation (|5.2|) to be 



2 r 2r 2 V ; 

where, as in the two-dimensional case, a 2 /2 is the strength of the 1/r 2 potential. The 
turning-point radii, determined from equation (|5.3|), are given by equations fl4.2|) and ( |1.3| ) 
with L 2 being the square of the angular momentum of the particle. As in the two-dimensional 
case, we choose the radius of the Poincare sphere to be R — 1/(2\E\). 

The calculation of the radial contribution to the action proceeds exactly as in two di- 
mensions. The result, similar to equation ( |4.5|) , is 



SUL 2 = (l + l) 2 h 2 ; E) = ir (^-) " - 2vr[(/ + \) 2 h 2 ± a 2 f 2 . (6.2) 

Setting the Maslov index o\ equal to 2, corresponding to two soft turnarounds in the radial 
motion (see Appendix A), we obtain from equation (|5.24|) the following condition for an 
energy eigenvalue for given I and m: 

/ 2 \ 1/2 

vr ( — J -2ti[{1 + \) 2 Ti 2 ±a 2 ] l/2 = 2vr^(n r + |), n r = 0,1,2,... (6.3) 

Since E is negative for the bound-state solutions we are considering, the energy eigenvalues 
for given / and m (—1 < m < I) are found to be 

E lmnr = - 2h 2 {nr + ! + [(/ + 1)2 ± a 2/ft2 ]1/2}2 > n r = 0,1,2,... (6.4) 

This expression gives the approximate semiclassical energy eigenvalues for the Coulomb plus 
1/r 2 potential. As expected, they are {21 + l)-fold degenerate because of the spherical sym- 
metry of the system. It should also be observed that the allowed values of I are constrained 
by the condition that {I + \)Ti < L max (E), with [L max (E)} 2 = 1/(2\E\) =f a 2 , as in the 
two-dimensional case. 

The result for the pure Coulomb potential is found by letting a — > in equation (p.4|). 
Putting n = ! + n r + lwe obtain, 

En = ~, n=l,2,... (6.5) 

independent of the sign of the a 2 term, as in two dimensions. This is the same as the 
familiar result found by solving the three-dimensional Schrodinger equation for the Coulomb 
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potential. It is clear from the definition of n that I < n. This condition, which also arises in 
solving the Schrodinger equation, leads to the degeneracy of the nth energy level being n 2 . 

It is remarkable that we have obtained the same result as the exact energy levels of the 
hydrogen atom, despite having made three approximations: (i) the semiclassical approx- 
imation embodied in Bogomolny's transfer operator; (ii) the asymptotic expansion ( |5.13| ) 
for Pz(cos7); and (iii) the evaluation of the integral in (|5.12| ) using the stationary phase 
approximation. A plausible explanation for this agreement will be given in the discussion 
at the end of the paper. 



B. The spherical harmonic oscillator plus 1/r 2 potential 

In this subsection we treat the isotropic harmonic oscillator plus a small 1/r 2 potential. 
We take the Hamiltonian corresponding to equation (5.2) to be 

tt Pr , 1 2 2 , L 2 ±a 2 

H = - + -ur +^^, (6.6) 

where, as in the two-dimensional case, u 2 describes the steepness of the harmonic oscillator 
potential, and a 2 /2 is the strength of the 1/r 2 potential, which may be attractive or repulsive. 
The radial turning points, determined from equation ( |5.3j ) for fixed E and L 2 , are given by 
equation ( |4.10| ), as in two dimensions. When L 2 has its maximum possible value, determined 
by 

[L max (E)] 2 = ^ T a 2 , (6.7) 

the radial kinetic energy is zero and the particle trajectory is confined to the sphere of radius 
R = E 1 ! 2 juj. We take the Poincare sphere to have this radius since all trajectories having 
L < L max (E) must repeatedly cross this surface. 

For given E and L 2 the radial part of the action can be calculated as in the two- 
dimensional case. The result is the same as equation ( f4.11|) . Furthermore, from the analysis 
in Appendix A, the Maslov index is o\ = 2. Thus, from equation ( |5.24|) , the condition for 
an energy eigenvalue is 

ttE 

Tih[(l + \f ±a 2 /h 2 } 112 = 2vr^(n r + |), n r = 0,1,2... (6.8) 

UJ 

Hence, the energy eigenvalue specified by I, m and n r is 

Eimn r = hu{2n r + [{I + \) 2 ± a 2 /h 2 }^ 2 + 1}, n r = 0,1,2,... (6.9) 

These eigenvalues are clearly (21 + l)-fold degenerate. Note that the permissible values of I 
are determined by the condition that (I + ~)h < L max (E), with L max (E) given by equation 

The energies of the pure isotropic harmonic oscillator in three dimensions are found by 
letting a — > 0. Defining n = 2n r + I, we obtain (independent of the sign of the a 2 term) 

E n = ?ku(n+ 1), n = 0,l,2,... (6.10) 
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The ground state energy is the zero-point energy associated with 3 freedoms. The multiplici- 
ties of the levels are determined by the (2/ + l)-fold degeneracy associated with m and by the 
number of distinct ways of obtaining a given value of n from the integer values of / and n r . 
The energy levels given by (|6.10|) and their degeneracies are the same as those obtained in 
the solution of the three-dimensional Schrodinger equation for the spherical oscillator (see, 



for example, Ref. [p5| , pp 166-168). 



C. Billiard in a spherical cavity 

A particle moving in zero potential inside a spherical cavity of radius R is the three- 
dimensional analogue of the circle billiard. As in the two-dimensional system, for given 
values of the energy E and the square of the total angular momentum L 2 , the inner turning 
point radius r_ is given by \L\ = {2E) l / 2 r_. We choose the PSS to be the sphere of radius R 
(or just slightly less than R, in order that the trajectories cross the PSS immediately after 
colliding with the spherical boundary). 

To obtain the EBK eigenvalues from equation Q5.24 ) we can use equation ( 4 . 1 8|) for 



the radial part of the action integral, which is the same for the two and three-dimensional 
systems. Furthermore, it was shown in Appendix A that the Maslov index is 07 = 3 for this 
system. Thus, from equations ( |5.24|) and (|4.18j ), the condition for an energy eigenvalue of 



the billiard in a spherical cavity is 

[k 2 R 2 -{1 + |) 2 ] 1/2 - {I + \) cos- 1 ^ + \)/{kR)\ = 7r(n r + |), 

1 = 0,1,2..., n r = 0,1,2,... (6.11) 

Here the allowed values of / at energy E are restricted by the condition (I + ~) < kR, 
where k = {2E) l l 2 /Ti. Equation (|6.11D could be solved numerically in a manner similar 



to equation ( |4.19| ) to obtain the EBK energy eigenvalues. (This would be equivalent to 
finding approximate values for the zeros of the spherical Bessel function j m (kR) when it is 
approximated by the leading term of the Debye asymptotic expansion.) 

VII. DISCUSSION 

We have shown how, with the help of the stationary phase approximation, Bogomolny's 
transfer operator leads to the EBK quantization rules for the energy eigenvalues of integrable 
systems having rotational symmetry, in both two and three dimensions. An important aspect 
of the theory was showing that the Maslov indices are correctly given by the simple rules of 
counting 1 for each soft turnaround and 2 for each hard-wall collision occurring during one 
complete cycle of the radial motion. 



In discussing the annulus billiard in Sec. IV E we drew attention to the fact that the EBK 



energy eigenvalues are appreciably different from those obtained using Bogomolny's transfer 
operator unmodified by making the stationary phase approximation or any other approx- 
imation. It is noteworthy, however, that among the eight different rotationally invariant 
systems to which the theory was applied, in five cases the EBK energy eigenvalues turned 
out to be the same as those obtained from an exact solution of the Schrodinger equation. 



In fact, for the hydrogen atom in three dimensions (Sec. [VI A| ), it was pointed out that the 
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correct energy levels were obtained despite having made three significant approximations in 
the theoretical development. How can one understand this surprising result? 

A similar situation arose in a recent paper |23[] concerning the application of Bogomolny's 
T-operator to a circular harmonic oscillator plus 1/r 2 potential. For that system it was 
possible to write down an exact transfer operator, which led to the exact energy eigenvalues 
without making any approximations. It was then shown, with the help of the Poisson 
summation formula, that the exact transfer operator could be written as an infinite sum of 
certain integrals, and that the leading term in this sum was the same as Bogomolny's transfer 
operator. Further, it was shown that improving the stationary phase approximation leads 
to corrections to the energy eigenvalues that involve higher powers of h than the leading 
term. In short, the semiclassical result agreed with the exact quantum result simply because 
corrections to the semiclassical approximation (which are nonzero) and to the stationary 
phase approximation (which are also nonzero) were not evaluated. Presumably, if corrections 
to all the approximations were evaluated systematically, they would cancel each other in all 
orders of h. 

It seems likely that we are dealing with a similar situation for the hydrogen atom in two 
and three dimensions, and for the harmonic oscillator with a singular magnetic flux line. 
For these systems it does not appear to be easy to write down an exact transfer operator, 
so that we are unable to carry out an analysis similar to that of Ref . [f23| . Nevertheless, if 
the EBK energy eigenvalue is regarded as the leading term of a semiclassical expansion in 
increasing powers of H, it is perhaps not so surprising that, in some cases, it agrees with the 
exact quantum result. 
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APPENDIX A: EQUIVALENCE OF a WITH THE MASLOV INDEX IN EBK 

QUANTIZATION 

In this Appendix we consider the index a m = fi m + v m introduced following equation 
( p,15| ) for a particle in two dimensions, and the index 07 = ^ + v\ defined in equation (|5.25|) 
for a particle in three dimensions. In both cases we consider a particle confined inside a 
finite region either by a smooth potential or a hard-wall boundary. From the definitions of 
the \i and v indices we show that a m and <j\ have the values one would compute from the 
simple rules for obtaining the Maslov index in EBK quantization. 

For two-dimensional systems, consider a particle moving in a smooth potential V(r). 
Figure 1 shows typical trajectories along with trajectories having slightly larger and slightly 
smaller angular momenta than the main trajectory. The trajectory with smaller angular 
momentum approaches more closely to the origin, but on the outward leg also ventures 
farther away from the origin. In the course of doing so, it must necessarily intersect the 
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trajectory with larger angular momentum. Thus, there is a focal point and /i m is thereby 
incremented by 1. There are then two possibilities. Before returning to the Poincare circle, 
the two orbits may intersect yet again, leading to \i m = 2, as shown in Fig. 1(a). If this 
happens one observes that the angle traversed is greater for the larger angular momentum, 
so that d 2 S/d , ~f 2 > and v m = 0. On the other hand, there may be no further intersection 
so that fi m = 1, as illustrated in Fig. 1(b). In this case one observes that d 2 S/d r y 2 < 
so that v m — 1. In either event we have a m = 2, which is the same as the value for the 
Maslov index in EBK quantization obtained by counting 1 for each of the turning points of 
the radial motion. 

We have performed numerical studies for two-dimensional systems having potentials of 
the form V(r) ~ ±r fc , where the plus sign is assumed for k > and the minus sign for k < 
(so that the potential is attractive). Our calculations show that for — 1 < k < 2 one finds 
fi m = 2 and v m = 0, while for —2 < k < — 1 and 2 < k < oo one finds fi m = 1 and v m = 1. 
The cases k — — 1 and k = 2, which are self- focusing, are marginal for the present analysis. 

The case of a hard wall may be approximated by taking the limit k — > oo. However, 
when considering both large k and infinitesimally close trajectories, one must be careful 
about the order in which the limits are taken. For any fixed k, no matter how large, we can 
take the two nearby angular momenta close enough to the main trajectory that we obtain 
a rn = 2, in conformity with the previous discussion. On the other hand, if we consider two 
fixed angular momenta, no matter how close, we can make k large enough that the two 
trajectories no longer intersect on the outer loop. This is the appropriate analysis for a 
disk with infinitely hard walls. In this case /i m = 2 from the collision with the hard wall 
(Dirichlet boundary conditions), while v m = 1, giving the result a m = 3. This agrees with 
the result for the EBK Maslov index obtained by counting 1 for the soft turnaround at the 
inner turning point and 2 for the hard turnaround at the disk boundary. The same analysis 
applies even if the motion is not force-free within the disk. For example, one could include a 
uniform magnetic field, or a flux line, or a harmonic potential out to the disk radius. As long 
as the trajectory bounces off the disk boundary, a m = 3, but if the particle avoids colliding 
with the boundary, the discussion in the previous paragraph will apply and a m = 2. 

Turning now to three-dimensional systems with a smooth potential, we recall that the 
Maslov index o\ was defined in equation ( |5.25| ) for the two different classes of trajectories. 
We now show that, for soft potentials, this leads to the result o\ = 2 in all cases. 

For given points Q and Q' on the Poincare sphere, a trajectory from Q to Q' lies in the 
plane containing these points and the origin. Focal points arising from variations within the 
plane of this trajectory yield the result [L\ + v\ = 2 from the preceding analysis. In three 
dimensions, however, there may exist an additional focal point lying on the straight line from 
the point Q to the origin. One can see that, provided the particle traverses an angle between 
7r and 2ir in moving along the trajectory for specified E and L 2 , it will intersect this line at a 
point on the opposite side of the origin from Q. Similar trajectories in all planes containing 
Q and the origin will intersect at this point, which is, therefore, a focal point. In this case, 
fii must be incremented by 1, giving \xi + vi = 3. But the angle traversed by the particle in 
this case is £j = 2n — ji, which means the trajectory belongs to the class £j = 2nji — with 
ji = 1. Hence, from equation (|5.25| ), the Maslov index is o\ = 2. Note that if the particle 
makes an additional circuit around the origin (corresponding to £j = 2nji — 7; with ji = 2) 
there will be two additional focal points on the line joining Q to the origin (on both sides of 
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the origin), but from (|5.25|) the Maslov index will still be 07 = 2. Clearly this generalizes to 
any number of complete circuits around the origin within the time of one cycle of the radial 
motion. 

If the angle traversed by the particle in moving along the trajectory specified by E and 
L 2 lies between and 7T, there is no focal point of the type described in the preceding 
paragraph. The angle traversed is £j = 7/ + 2nji with ji = 0, and equation (|5.25|) gives 

01 = 2. If there are additional circuits around the origin corresponding to £j = 7; + 2nji, 
with ji = 1,2,..., one can easily see that /ii must be incremented by 2 for each circuit, but 
from equation ( p.25| ), the Maslov index remains oi = 2. Thus, we have shown that 07 = 2 for 
all possible trajectories in a soft potential. This result is the same as would be obtained from 
the simple EBK rule of counting 1 for each radial turning point of the effective potential 
V(r) + L 2 /(2r 2 ). 

If the outer radial turning point is replaced by a hard-wall collision with the boundary 
(Dirichlet boundary condition on the wave function), the above analysis still holds up to 
the point of colliding with the boundary. As for the two-dimensional systems, the hard-wall 
collision requires incrementing in by 1, giving the result o\ = 3. This is the EBK result from 
the simple rule of counting 1 for the soft turnaround at the inner radial turning point plus 

2 for the collision with the boundary. 



APPENDIX B: THE DETERMINANT OF SECOND DERIVATIVES OF S 



In this appendix we evaluate the determinant of second derivatives of the action, which 
constitutes the amplitude of the three-dimensional transfer operator: 



det- 



d 2 S 



d 2 S 

86 86' 
8 2 S 

86 86' 



8 2 S 

86 86 
8 2 S 



(Bl) 



Here S = S(£±;E), which, through £^ , depends on the angle 7 introduced just after 
equation Q5.4|). Since 7 is the angle subtended at the origin by the points Q and fl' on the 
Poincare sphere, we have 



cos 7 



sin 6 sin 6' cos(6 — 6') + cos 6 cos 6' 



< 7 < ?r. 



(B2) 



Our objective is to express the determinant in terms of 7. 
First, let us write 



8S_ 
~86 



as 8j_ 
97 86 



8 2 S 
86 86' 



(PS_ 8j_ 8j_ 8S_ (9 2 7 
St 1 86 86 1 + ~8^ 86 86' ' 



(B3) 



It will greatly simplify the calculation to evaluate the second derivatives assuming that 
6 = 6' = 7r/2. This means that the polar axis of the spherical polar coordinate system 
is chosen to be perpendicular to the plane containing the origin and the arbitrarily chosen 
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points (R, Q) and (R, Q') on the PSS. Because of the spherical symmetry, the result obtained 
will be valid for any orientation of the axes of the spherical polar coordinate system. When 
= 6' = it/2, one finds that d^/dO = fry/dP = and d 2 ^/(d9d9') = -l/sin 7 . Hence, 



d 2 S \ 1 dS 



\d6d8' J e=e , =7T/2 sin 7 d-f 

The other second derivatives can be evaluated in the same way. We obtain 

{ d 2 S \ ( d 2 S \ 



d9d<P') e=el=w/2 \d<pd6> ) e=e , =7r/2 
( d 2 S \ d 2 S 







\d4>d4>'J d=dl=n/2 <9 7 2 
Hence, the determinant in (|B1|) has the value given in equation (|5.5|) . 



(B4) 



(B5) 
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FIGURES 



Examples of trajectories in two dimensions which start from an arbitrary point on the 
Poincare circle and return (in the same sense) to the Poincare circle at the point 0'. In each 
case there is a main trajectory together with trajectories having slightly larger and slightly 
smaller angular momenta. Fig. 1(a) was computed for the potential V(r) ~ r k with k = 0.5. 
Fig. 1(b) was computed for V(r) ~ r k with k = 3.5. 
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